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Intrigued by a recent Belle result for a large direct CP asymmetry in B° — > 
D + D~, we study the effects of a b — > wwc/ quark transition by combining 
the asymmetry information with rates and asymmetries in isospin-related 
decays. Arguing for a hierarchy among several contributions to these decays, 
including an exchange amplitude which we estimate, we present tests for 
factorization of the leading terms, and obtain an upper bound on the ratio 
of b — > uud and b — > ccd amplitudes. We prove an approximate AI =1/2 
amplitude relation for B — ► DD, and an approximate equality between CP 
asymmetries in B° — > D + D~ and B + — > D + D°. Violations of these relations 
by Belle measurements, at 1.8a and 3.6cr respectively, if confirmed, would 
indicate a possible New Physics contribution in b — > uud. Applying flavor 
SU(3), we extend this study to a total of ten processes including AS = 
decays involving final D s and initial B s mesons, and AS = 1 decays of B and 
B s mesons into pairs of charmed pseudoscalar mesons. The decays B s — > DD 
provide useful information about a small exchange amplitude, responsible 
for a decay rate difference between B + — > D + D° and B° — > D + D~ . A 
method for determining the weak phase 7, based on CP asymmetries in 
B°(t) -> D + D~ and the decay rate for 5 S -> £>+£>7 or B+W -»■ D+D°(D~), 
is shown to involve high sensitivity to SU(3) breaking. 

PACS codes: 12.15.Hh, 12.15.Ji, 13.25.Hw, 14.40.Nd 



1 Introduction 

Accurate measurements of the weak phase (3 = a,Tg(y t * b Vtd/V* b V c d), sin 2/3 = 0.680 ± 
0.025, cos 2(3 > [1], have provided a precision test for the Cabibbo-Kobayashi-Maskawa 
[2,3] framework and for the Kobayashi-Maskawa mechansim of CP violation. The accu- 
racy of this test relies on the pure dominance by a single weak phase of a few b — ► ccs 

1 On sabbatical leave from the Physics Department, Technion, Haifa 32000, Israel. 
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Decay mode 




BaBar 


r>elle 


Average 


B° — > D + D~ 


B 


28±04±0 5 


1 97 ± 20 ± 20 


2 11 ± 31 




A C p 


-0.11 ±0.22 ±0.07 


0.91 ±0.23 ±0.06 


0.37±0.17 




S 


-0.54 ±0.34 ±0.06 


-1.13 ±0.37 ±0.09 


-0.75 ±0.26 


B + -> 


B 


3.8 ±0.6 ±0.5 


3.85 ±0.31 ±0.38 


3.84 ±0.42 




A C p 


-0.13 ±0.14 ±0.02 


0.00 ±0.08 ±0.02 


-0.03 ±0.07 




B 


< 0.6 (90% c. 1.) 


< 0.43 (90% c. 1.) 


< 0.43 (90% c. 1.) 




B 


< 1.0 (90% c. 1.) 


< 0.36 (90% c. 1.) 


< 0.36 (90% c. 1.) 



Table I: Charge-averaged branching ratios B in units of 10 4 and CP asymmetries Acp, S 
in B — > DD, from Refs. [10-13]. Also included are upper limits on B(B° — > DfDj) [14, 
15]- 



processes including 5° — > J/ipKs [4,5]. This implies a mixing-induced asymmetry, 
S* = sin 2/5, and a vanishingly small direct CP asymmetry, as confirmed experimentally, 
A CP = -0.012 ±0.020 [1]. 

The decay B° — > D + D~ is dominated by b — > ccrf, but involves a smaller non- 
negligible amplitude from 6 — > wwd carrying a different weak phase. The second ampli- 
tude introduces hadronic uncertainties in predictions for the asymmetries S and Aqp 
in this process. Early model-independent estimates of the ratio of the two amplitudes 
contributing to this process vary from a few percent to upper bounds of about 0.2 [4,6] or 
0.3 [7]. A more recent model-dependent calculation finds 0.03 [8]. Values larger than 0.3 
may be obtained in extensions of the Standard Model [9] . The two asymmetries depend 
also on the strong phase difference between the b — > ccd and b — > uud amplitudes [4] . 

Asymmetries in B° — > D + D~, measured by the BaBar and Belle collaborations, 
are quoted in the upper part of Table I. The table also includes branching ratios for 
B° -> D+D-,B + -> D + D°,B° -> D°D°, a direct CP asymmetry measured for B + -> 
D+D° [10-13], and upper limits on B(B° -> D+D;) measured by BaBar [14] and 
Belle [15]. The BaBar asymmetries in B° — > D + D~ [10] are consistent with Aqp = 
0, S = — sin 2(3, showing no evidence for a b — > w-u<i term in the decay amplitude. In 
contrast, the Belle asymmetry measurements [11], which fluctuate outside the physical 
region, A 2 CP + S 2 < 1, deviate substantially from the above nominal values, indicating a 
sizable b — > uud amplitude. The direct asymmetry Aqp measured by Belle is nonzero at 
a level higher than 3o\ Its central value indicates the possibility of a second amplitude 
larger than permitted in the Cabibbo-Kobayashi-Maskawa (CKM) framework. 

A major goal of this paper, largely intrigued by the Belle results, is to study carefully 
the dynamics and CKM structure of the B° — > D + D~ decay amplitude and of decay 
amplitudes for the two isospin-related processes, B + — > D + D° and B° — > D°D°. In 
references [8] and [16] these processes have been stated to originate in a AI = 1/2 
effective Hamiltonian implying an isospin triangle relation among the three amplitudes. 
It will be shown that, while AI = 1/2 is not a property of the effective Hamiltonian, 
an approximate AI = 1/2 rule is expected to hold for the three decay amplitudes and 
should be tested experimentally. Applying flavor SU(3) to the above processes, we will 
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extend our study to include strangeness-conserving decays involving final D s and initial 
B s mesons, and strangeness-changing decays of B and B s mesons into pairs of charmed 
pseudoscalar mesons. 

In Section 2 we study the asymmetries measured by BaBar and Belle in B° — > D + D~ 
in terms of two parameters, the ratio r of b — > uud and b — > ccd amplitudes and their 
relative strong phase 5. Section 3 introduces expressions for the amplitudes of the three 
processes B° — > D + D~,B° — > D°D° and B + — > D + D° in terms of isospin amplitudes. 
We identify circumstances under which an isospin triangle relation between these ampli- 
tudes can be violated by a (small) AI = 3/2 contribution. Section 4 studies B — > DD 
decays and two other SU(3) related AS = decays of B° and B s in terms of graphical 
contributions, while Section 5 extends this study to corresponding strangeness changing 
decays of B°, B + and B s . Section 6 discusses a hierarchy among graphical amplitudes, 
presenting tests of factorization for the dominant terms. In Section 7 we discuss briefly 
consequences of this hierarchy on a theoretical upper limit on r, illuminating an inconsis- 
tency between the CP asymmetries measured by Belle in B° — > D + D~ and B + — > D + D°. 
Section 8 discusses a way for determining the weak phase 7 by combining information 
from asymmetries in B° — > D + D~ and decay rates of corresponding AS = 1 decays 
to charm-anticharm, while Section 9 concludes. An Appendix provides a dictionary be- 
tween graphical amplitudes and SU(3) reduced matrix elements of four-quark operators 
appearing in the effective Hamiltonian. 



2 Ratio of b — > uud and b — ► ccd terms in B° — > D + D 



We start our discussion by translating the B° — > D + D~ asymmetries, measured sepa- 
rately by BaBar and Belle, into values of the ratio r of b — > mmJ and b — > ccd amplitudes 
and the relative strong phase 5 between these amplitudes. Denoting 



A(B° -> D+ZT) = A c + A u e i(<5+7) = A c [l + r e i(<5+7) ] , (r = A u /A c ) , 




(1) 



one has [4] 



Acp(D + D-) 



S(D + D-) 



2Im(A D+D -) 
1 + |A D+D -| 2 

\^D+D- 1 2 ~ 1 

|A D+D -| 2 + 1 



sin 2/? + 2r cos 5 sin(2/3 + 7) + r 2 sin 2(/3 + 7) 
1 + 2r cos 5 cos 7 + r 2 
2r sin <5 sin 7 



(2) 



1 + 2r cos 5 cos 7 + r 2 



Keeping only linear terms in r, 



S(D + D-) 
A CP (D + D-) 



~ — sin 2/3 — 2r cos 2(3 cos 5 sin 7 , 
~ 2r sin 5 sin 7 , 



(3) 
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Figure 1: \ 2 plots in the (r cos 5, r sin 5) plane assuming (3 = 21.5°, 7 = 68°. The red, blue 
and green curves show constraints following from B° — > D + D~ asymmetries measured 
by BaBar, Belle and their averages. In each case the most inside, intermediate and most 
outside curves represent bounds at 68%, 90% and 95% confidence levels. Red and green 
points describe solutions corresponding to central values of the BaBar asymmetries and 
the averaged asymmetries. 



Consider the measured asymmetries and the current values of (3 and 7, (3 = (21.5 ± 
1.0)° [1], 7 = (67.6l|;|)° [17] (see also [18,19]). Using this information, the approximation 
(J3j) and (j4lh or the precise expressions (j2J), determine r and 5. The resulting errors in r 
and S are dominated by the errors in the measured asymmetries. Taking central values 
for the asymmetries and values (3 = 21.5°, 7 = 68°, Eq. (j3J) implies central values around 
r = 0.1 (BaBar) and r = 0.6 (Belle). In both cases the error in r is about 0.2. The 
central value of r for Belle, for which the linear approximations ([3]) involve non-negligible 
quadratic corrections, should be considered with care because this value of r is based on 
non-physical values of the asymmetries obeying A^ P + S 2 > 1. 

In order to study the implications of the asymmetries on the pair of parameters 
(r, 5), we have performed a two dimensional \ 2 analysis for these two parameters using 
the asymmetry measurements and assuming (3 = 21.5°, 7 = 68°. In Fig. [TJ we plot the 
resulting contours in the (r cos S, r sin 5) plane for x 2 — 2.30,4.61,5.99, corresponding 
to 68%, 90%, 95% confidence levels. The red, blue and green curves show constraints 
following from the asymmetries measured by BaBar, Belle and their averages. In each 
case the innermost, intermediate and outermost curves describe bounds at 68%, 90% and 
95% confidence levels. The red and green points are solutions corresponding to central 
values of the BaBar asymmetries and the averaged asymmetries. We do not show a 



implies 




(4) 



2 sin 7 
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central Belle point because the central values of the Belle asymmetries lie outside the 
physical region. 

The Belle two-parameter boundary curve for 90% confidence level contains a point 
with closest distance to the origin, r = 0.29,5 = 82°. The point on this curve with 
smallest 5 has 5 = 37°, r = 0.87. (Note that the values r = 0.29 and 5 = 37° are lower 
than the 90% confidence level lower limits on these separate single variables.) Thus, 
the Belle data alone would provide evidence for a sizable r and for a large strong phase 
difference 5. We note, however, that the BaBar and Belle regions of 90% confidence level 
do not overlap. In order to draw firm conclusions about r and 5 one should therefore 
wait for better agreement between the asymmetries measured by the two collaborations. 



3 Isospin amplitudes in B — > DD 

The low energy effective Hamiltonian governing B° — > D + D~ , B° — > D°D° and B + — > 
D + D° involves two CKM factors V* b V cd and V* b V ud , both of order A 3 (A = \V US \ = 
0.2258 ± 0.0010 [17]). Each of these factors multiplies a combination of four quark 
operators with coefficients given by calculable Wilson coefficients Cj [20], 



n-f - Gp 



2 2 10 

v: b v cd £ ca c + v: b v ud £ dot + {v: b v cd + v: b v ud) £ c k o k 

i=l i=l fc=3 



(5) 



The current-current operators 0\ and 0\ (i = 1,2) have flavor dependence (bc)(cd) 
and (bu)(ud), respectively. Thus, while the first pair of operators are pure AI = |, 
the second pair involves both AI = | and A/ = |. The QCD penguin operators O k 
{k — 3 — 6) with flavor structure are pure AI = 1/2, while the electroweak penguin 
operators O k ~ (pd)Y,q e qisi ( l) (k — 7 — 10, g = u,d,s,c), which depend on the quark 
charges e g , involve both AI = 1/2 and AI = 3/2. Thus, in contrast to statements 
made in Refs. [8] and [16], the effective Hamiltonian underlying B — > DD decays is not 
pure AI = 1/2. It contains an additional AI = 3/2 component, also when neglecting 
electroweak penguin contributions which have very small Wilson coefficients [20]. 

The final DD states consist of / = and 1=1. This implies a total of three 
isospin amplitudes, A i,A 1 i and A^, where the two subscripts denote I(DD) and 
AI, respectively. Neglecting very small electroweak penguin contributions, the AI = 3/2 
amplitude A ± ^ occurs in association with a CKM factor V* b V ud but not with V* b V cd . 

The three physical B — > DD decay amplitudes can be written in terms of the three 
isospin amplitudes, 

A + - = A(B° — > D + D~) = - A, i + - A i + - A s , 
v 2 2 2 2 2 1 '2 ' 

A +0 = A(fi + ^/J + /3°) = A a -iA lif . (6) 
These relations can be inverted, 



i oo 
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A ll = ~(A + - + A 00 + 2A +0 ) , 

' 2 3 

A,s = 2 -(A+~ + A™ - A+») . (7) 

While AI = 1/2 is not a property of the low energy effective Hamiltonian, we will 
argue below that \A 1 3\ <C \A i 1, \A l i 1 is a reasonable assumption which should be 
tested experimentally. Neglecting the AI = 3/2 amplitude, one obtains an approximate 
triangle relation [8,16], 

A+- + A 00 = A +0 . (8) 

A potential proof for a nonzero A J = 3/2 amplitude would be a violation of ([8]). This 
happens when the amplitude triangle does not close, for instance when + \ A 00 \ < 

\A»\. 

In order to illustrate such a possibility consider the Belle measurements for B — > DD. 
We define 



\A + -\ = 10 2 JB(B° -> D+D-) [1 - A CP (B° -> D+D-)} 



|A +0 | = 10 2 \jB(B+ -> D + D°) [1 - A CP (B+ -> D+Z)°)J (r /r+) , 

|A 00 | = 10 2 -i/i3(5 -»• D°D°) \l - A CP {B° ->■ D°Z) )j . (9) 



Using the Belle values for CP-averaged branching ratios and CP asymmetries given in 
Table U and a ratio of B + and 5° lifetimes [1] r + /r = 1.071 ± 0.009, we find 

\A + -\ = 0.42 ±0.56 , |A +0 | = 1.90 ±0.14 , \A 00 \ < 0.57 (la) . (10) 

The small magnitude of A + ~ follows from the large positive CP asymmetry measured 
by Belle, implying observing mostly B° decays into D + D~ with only a few B° decays 
into this final state. The upper bound on |y4 00 | is obtained from a la upper limit on 
B(B° — > D°D°). We have assumed that the CP asymmetry in B° — > D°D° is not large 
and negative (in Section 7 we will argue for a vanishingly small asymmetry), and we 
have neglected possible correlations between errors in branching ratio and asymmetry 
measurements in the other two modes. We note that the triangle ([8]) does not close for 
central values of the measured amplitudes fllOp . and allowing for deviations from these 
values up to 1.8cr . 

In contrast, the triangle relation holds well when using the central values of the 
BaBar measurements. A critical test for the closure of the amplitude triangle requires a 
reduction in errors and a better agreement between BaBar and Belle. 



4 AS = decays into DD using graphical amplitudes 

Useful expressions for amplitudes in B — > DD using flavor SU(3), which can be gen- 
eralized to B° — > D+D~ and B s — > D + D~, are in terms of graphical contributions 
representing the flow of isospin and flavor SU(3) in these decays [21]. This includes 
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a (color-favored) tree amplitude T involving V* b V c d, penguin and penguin- annihilation 
amplitudes P and PA with u, c and t-quark loops, exchange amplitudes E c and E u 
involving V* b V c d and V* b V u d, and an annihilation amplitude A u involving V* b V u d- The 
amplitude E c is associated with popping a pair of uu or dd out of the vacuum, while 
E u and A u involve cc popping. At this point we neglect very small electroweak penguin 
contributions to which we return in the next section. 

The graphical amplitudes have well-defined isospin properties. The two graphical 
amplitudes E u and A u involve both AI = 1/2 and AI = 3/2, while all other amplitudes 
are pure AI = 1/2 by construction. We will show below that the AI = 3/2 amplitude 
consists solely of the combination E u + A u . This combination is thus responsible for a 
potential violation of the amplitude triangle relation (jSJ). 

We denote 

p = v: b v ud p u + v: b v cdPc + v*v tdPt = v: b v cdPct + v* b v udPut , ( Pij = Pi - Pj ) , (11) 

absorbing the first term in the tree amplitude by defining, 

v: b v cd t c = T+v: b v cdPct . (12) 

Similarly, writing 

PA = V* b V ud P a u + V* b V cd P a c + V* b V tdP a t = V* b V cd P a ct + V* b V ud P a ut , (pay = P a { - P a j ) , 

(13) 

the first term is absorbed in E c , 

V; b V cd e c = E c + V: b V cd P a ct . (14) 
Other terms involving the CKM factor V* b V u d are 

E u = V* b V ud e u , A u = V* b V ud a u . (15) 
Using these definitions with a shorthand notation, Pu = Pu t,pa u = pa>ut, w e find: 



a. 


A(B° - 


> D + D~) 


= v; b v cd 


(t c + e c ) + V* b V ud ( Pu + pa u 


b. 


A(B° - 


-> D°D°) 


= v; b v cd 


(-e c ) + KTfeKd (-pa« - e u ) 


c. 


A(B + - 


-> D+D°) 


= v; b v cd 


(t c ) + KfeKd (p« + o u ) , 


d. 


A(B° - 


♦ DtD-) 


= v: b v cd 


( e c) + V* b V ud ( P a u ) , 


e. 


A(B S - 


> d+d;) 


= VSYa 


(t c ) + V* b V ud (p u ) . 



(16) 

The minus sign in the amplitude involving a -D follows from our convention, D° = 
-cu [21]. 

The five physical amplitudes depend on the four combinations, V* b V cd t c + V* b V udPu , 
V* b V cd e c + V* b V ud P a u , V* b V ud e u , V* b V ud a u . This equals the number of independent SU(3) 
reduced matrix elements describing AS" = B and B s decays to pairs of charmed 
pseudoscalar mesons, (1||3||3), (8||3||3), (8||6||3), (8||15||3) (see Appendix A). Consequently 
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\Pu +pa u \ 


\v: b v cd \ 







the five decay amplitudes are not mutually independent. They obey one triangle relation 
in the SU(3) symmetry limit, 

A{B° -> D+D~) + A(5 S -> £+£;) = -> D+ZT) . (17) 

The parameters r and 5 studied in Section [2] can be expressed in terms of graphical 
amplitudes, 

/ Pu+pa u 
arg — — 

V t c + e c 

and the isospin amplitudes defined in section [3] are given by 

A) ± = v *b v cd (t c + 2e c ) + V* b V ud (p u + 2pa u + e u ) , 
A,i = V* b V cd t c + \v* h V ud (3 Pu -e u + 2a u ) , 

\l = -\v: b V ud (e u + a u ) . (19) 

The last relation confirms our statement above that the AI = 3/2 amplitude involves 
only the combination E u + A u . 

5 AS = 1 B and B s decays into charm- ant icharm 

The parametrization (fTSj) of AS = decays into charm- ant icharm in terms of flavor 
SU(3) graphical amplitudes can be extended to AS = 1 CKM-favored decays of B and 
B s mesons, which are governed by b — > ccs. U-spin reflection symmetry d <-> s implies 
expressions similar to Eqs. (fTE|) for corresponding U-spin related amplitudes [22], in 
which one replaces V* b V cd by V* b V cs and V* b V ud by V^&V^g. Since U-spin transforms 
B° ^ B s , D ± Df while keeping B + , D° and D° invariant, one has: 

a. A{B S -> D+D7) = V c ;Ks (t c + e c ) + V* b V us (p u + pa u ) , 

b. A(B S -> D°D°) = V; b V cs (-e c ) + V£K. (-pa„ - e u ) , 

c. A(5+ -> D+£>°) = V^V^ (* c ) + KbK, (p„ + a tt ) , 

d. A(5 S - D+D") = V; b V cs (e c ) + V: b V us (pa u ) , 

e. A(B° — > D+D~~) = V: b V cs (t c ) + V: b V us ( Pu ) . (20) 

We note that while both V^V^ and V* b V ud in ([To]) are of order A 3 , the factors V c * fe V cs 
and V^Vus in f[2"0l are order A 2 and A 4 respectively. More precisely, one has 



KTfeKs A 



V* b V ud 1 - A 2 /2 ' 



(21) 



This implies that in the U-spin symmetry limit CP rate differences in the five U-spin pairs 
of corresponding AS = and AS = 1 decays, such as B° — > D + D~ and -B s — > D+D~ , 
have equal magnitudes but opposite signs [22], 

r(B s -> d+z>7) - r(fi s -> d+i>7) = -[r(B° -> D+zr) - r(#° -> z>+zr)] . (22) 



8 



Decay mode: 


B° -»■ D+D~ 


B + -> DfD [> 




Branching ratio: 


79 ±7 [15,24,25] 


109 ±18 [24,25] 


«CS [26] 



Table II: Charged averaged branching ratios in units of 10 4 for AS" = 1 B and B s 
decays into charm-anticharm. 

Since the rates of AS = 1 decays are about A~ 2 larger than those of AS = decays, the 
CP asymmetries in the former are expected to be about A 2 smaller than in the latter. 
The case of Eq. ( 122]) has been discussed in Ref. [23]. Similar CP rate difference relations 
hold between the other four pairs of processes. An SU(3) triangle relation analogous to 
(ITT)) holds for AS = 1 decays, 

A(B S -> D + D~) + A(B° -> D+D~) = A(B S -> D+Dj) . (23) 

So far we have neglected electroweak penguin contributions in both AS = and 
AS* = 1 decays. This is justifiable in the first case by very small Wilson coefficients 
associated with electroweak penguin operators [20] multiplying a CKM factor V* b V t d of 
the usual order A 3 . However, electroweak penguin terms in strangeness-changing decays 
involve a CKM factor V t * b V ts ~ 0(X 2 ), much larger than terms involving V* b V us ~ C(A 4 ) 
which are kept in (T2"0"|) . Thus, for consistency one must keep also electroweak penguin 
terms in AS" = 1 decays. 

In general, there are four types of diagrams describing electroweak penguin (EWP) 
contributions, corresponding to the above-mentioned four SU(3) reduced matrix ele- 
ments. The four diagrams may be associated with a color-suppressed EWP amplitude 
P EW , an EWP-exchange amplitude PE EW (c) associated with cc popping, and two EWP- 
annihilation amplitudes, PA EW and PA EW (c), associated with uu,dd,ss popping and 
cc popping, respectively. We neglect the two EWP amplitudes involving cc popping for a 
reason discussed in Section [6j Using V t * b V t d(s) = —V* b V c d( s ) — KTfeKtd(s), the remaining two 
EWP amplitudes, P EW and PA E w, may be absorbed in the following way in definitions 
of four amplitudes occurring in (|T6|) and (I2"0"j) without changing these ten equations: 

2 2 2 2 

t c--^P EW ^ t c e c --PA EW -> e c ,p u --P EW -> Pu , pa u --PA EW ^ pa u . (2A) 

We conclude this section by quoting in Table HT1 branching ratios for three of the five 
processes occurring in Eqs. (120]) . including a very recent measurement of B(B S — > D+D~) 
by the CDF Collaboration [26]. 

6 Expected hierarchy among graphical amplitudes 

6.1 AS = 1 decays 

6.1.1 Ratio of two CKM factors 

Consider first the five AS = 1 amplitudes ( |20l) for B and B s decays into pairs of 
charmed pseudoscalar mesons. Each amplitude involves a dominant term with a CKM 
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factor V* b V cs and a much smaller term involving V* b V us . The CKM suppression of the 
second amplitude, which is often being neglected, is [17]: 



ub y us\ 



\VlVc 



(0.40 ± 0.05)A 2 /(1 - A 2 ) = 0.021 ± 0.003 . (25) 



cb v cs\ 



6.1.2 The amplitudes t c and e c 

The large CKM factor V* b V cs multiplies two amplitudes, a dominant term t c and a smaller 
exchange contribution e c accompanied by an EWP-annihilation contribution PAew [see 
f l24|) ]. The latter two terms are suppressed by Aqcd/^-b as they involve the interaction 
of a light spectator quark. Thus, we expect the two decay modes B s — ► D + D~ and 
B s — ► D°D° governed by e c to have branching ratios much smaller [27] than those quoted 
in Table U for B° -> D+D-,B + -»■ D+Z)° and fi s -> £>+£>J which are dominated by 
£ c . The corresponding small ratios of branching ratios, e.g. B(B S — ► D + D~)/B(B° — > 
D+D~), would determine |e c /t c | 2 . 

A reasonable although not precise estimate for |e c /t c | may proceed as follows. Con- 
sider the two processes B° — > D~n + and B° — > D~K + , both of which originate in 
the quark subprocess b — > cud. While the first decay is governed by a color-favored 
tree amplitude with a small exchange contribution, the second one obtains only a con- 
tribution from an exchange amplitude. Drawing a parallel between these amplitudes 
and the corresponding amplitudes in b — > ccs, and using [25] B(B° — > D~n + ) = 
(26.8 ± 1.3) x 10" 4 , £(5° -> D-K+) = (2.8 ± 0.5) x 10" 5 , we estimate 



|t c + e c | ' " ^ i3(5° 
In this crude approximation this would imply 

le, ' 



0.102 ±0.009. (26) 



0.093 ± 0.008 < y^r < 0.114 ±0.010 . (27) 

\tc\ 

We note three corrections which may affect this estimate: 

1. The amplitude e c includes by definition a term PAew (both terms require an in- 
teraction of the spectator quark), while no such term contributes to B° — ► D~K + . 

2. The exchange amplitude in B° — ► D~K + involves popping in a cu system, 
whereas e c is described by uu or dd popping in a cc system. 

3. B° — > D~7i + is dominated by a purely a color-favored tree amplitude, while t c 
involves also a smaller penguin term p ct . See discussion below. 

These three differences are expected to affect (|27|) by a factor which is hard to estimate. 

In our discussion below we will make a conservative assumption based on measure- 
ments for AS = decays, 



B(B° -> D°D°) t + . . 

f < 0.3 . (28) 



t c | \ i3(fi+ -> D + D°) T 
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This upper bound is obtained from the branching ratios quoted in Table H for B° — > 
D°D° and B + — > D + D°, for which la upper and lower limits are used. The two 
processes B° — > D°D° and B + — > D + D° are dominated by e c and t c , respectively, while 
smaller contributions involving V* b V u d have been neglected. (See discussion below.) 
It would be useful to compare the bound fl28|) and the crude estimate fl27|) with direct 
measurements of \e c /t c \ in ratios of branching ratios for AS* = 1 decays including B(B S — > 
D + D-)/B(B + -> D+D°) and B(B S -> D°D°)/B(B -> D+D~). Using the averaged 
measured B(B + — > D + D°) in Table I to normalize |t c | 2 , Eq. (1271) would imply — > 
DfDj) = (4. Ol};^) x 10 -6 , about an order of magnitude smaller than the current upper 
limit on this branching ratio (see Table I.) 



6.1.3 The ratio |V^V cd p ct /T| in t c 

The amplitude t c consists of a combination of a genuine color-favored tree amplitude and 
a smaller loop-suppressed penguin amplitude, p ct , with t and c quarks in the loop [see 
Eq. (j!2p ]. It is difficult to obtain a good estimate for the ratio of these two amplitudes, 
the sum of which contributes to both AS* = and AS = 1 decays. A QCD loop 
factor [a s (m&)/127r] ln(m 2 /m 2 ) characterizing the suppression of V* b V c dp c t relative to T 
(or a typical Wilson coefficient for penguin operators [20]) is of order five percent. A 
dynamical enhancement by a factor of four to six relative to a QCD loop factor has been 
measured for the penguin-to-tree ratio in B° — > 7r + 7r~ [28,29]. We will permit a similar 
enhancement in B° — > D + D~ , thus allowing | V^V^Pctl to be at most as large as 0.3|T|, 



< 0.3 . (29) 



6.1.4 Factorization of (V cs /V cd )T in B° -> D+D~ and fi+ -> DfD° 

The tree amplitude is expected to factorize within a reasonable approximation into 
a product of a _B ^ D form factor and the D s meson decay constant. While this 
approximation cannot be justified by the heavy h quark limit of QCD (which can only 
be applied when a B meson decays into two energetic mesons [30,31]), it holds to leading 
order in 1/N C in the large N c limit [32]. Early factorization tests of this kind, implicitly 
neglecting a p ct contribution, have been proposed and studied in Ref. [33,34]. 

We will now update a factorization test for (V cs /V c d)T in B° — > DfD~ and B + — > 
DfD° by relating the branching ratios for these two processes given in Table [III to the 
above quoted branching ratio for B° — ► D~tt + . The latter has been accounted rather 
well by factorization [31], up to a small contribution from an exchange amplitude. The 
dominant contributions of the isosinglet amplitudes (V cs /V c d)T and V* b V cs p ct to the decay 
rates of B + — > D+D° and B° — > D^D~ are expected to be equal in the isospin symmetry 
limit. Thus we will use the weighted average of these two branching ratios, correcting 
B(B + -> D+D°) by the lifetime ratio r /r + : 

B(B -»■ £>+£>) = (82.4 ± 6.5) x 10~ 4 . (30) 
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This implies the following ratio of measured branching ratios, 



Bf^OiD) , 3 _ 07±a2g _ (31) 



B(B° -> L)-tt+ 

Assuming factorization for the processes in the numerator and denominator and 
taking V cs /V u d = 1, one would expect this ratio to be given by [34] 

B(B — > DfD) _ f 2 D3 F^D s )[(l + VC^ 2 -Cn s } 2 PD a 



B(B°^D-n+) ft F$M [(l + v^) 2 -^] 2 Ptt 

where 



(32) 



Hi + Hiz™! (33) 



Here /o s , and pn 3 , P-k are corresponding decay constants and momenta in the B 
meson rest frame, while Fy is the B —> D vector form factor. Taking a linear form 
factor [35], F v (u) = F v (l)[l - (0.69 ± 0.14) (w - 1)], and using [36] f w = 130.4 ± 0.2 
MeV, fo s = 273 ± 10 MeV, with meson masses and momenta in the B rest frame listed 
in Ref. [25], one finds 

B(B — > D+D) _ 1AM7 , 



Comparing the factorization calculation with the experimental ratio (13~TT) we note 
that the factorization result is on the high side, showing a discrepancy of 2.1<r relative 
to experiment. We do not expect a very good agreement here because of two corrections, 
each of which may be about 30% in amplitude: 

1. A penguin term V* b V cs p ct contributing to the numerator but not to the denom- 
inator. As mentioned, this term could be as large as 0.3T and could interfere 
destructively with T, leading to a ratio smaller than (I34p by up to a factor of two. 

2. Nonfactorizable 1/N C corrections to T contributing to both numerator and denom- 
inator. These terms are also expected to lead to corrections around 30% in the 
amplitude. 



6.2 AS = decays 

6.2.1 Ratio of two CKM factors 

We now turn to the AS = decay amplitudes given in Eqs. ffTBI) . These amplitudes 
involve the two CKM factors V* b V C( i and V* b V u d which are of comparable order A 3 , with 
a ratio [17] 

= 0.40 ± 0.05 . (35) 

\v; b v cd \ 1 ; 

In the U-spin symmetry limit these CKM factors multiply the same hadronic amplitudes 
occurring in AS = 1 decays. The somewhat larger CKM factor V* b V c d multiplies a 
dominant term t c and an exchange contribution e c which is expected to be much smaller. 
This leads to a large suppression of B(B° — > D°D°) and B(B° — > D+D~) relative to 
B(B° -> D+D-), B(B+ -> D+D°) and B(B S -> D+D;) [27, 37]. 
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6.2.2 Factorization of T in B° -> D + D~ and 5+ -> L> + L>° 

The hadronic amplitude V* b V cd t c consists of a dominant term T, which is factorizable in 
the large N c limit, and a sub-dominant term V^V^Pct, which may be as large as about 
0.3T. Using notations as above, factorization of T implies that the contribution of this 
amplitude to the rates for B° — ► D + D~ and 5 + — > D + D° are equal and are given by 



B(B° -> D + D-) T B(B + -> D+D 



OA 



o.i3(r;;;-. <.■>>(» 



A 2 /£ F* [(1 + ^ - Cd ]2 Pd 



" i-a 2 /2^k) [(i + VG) 2 - C.] 2 p. 

We have used a value [36] fo = 205.8 ± 8.9 MeV. Using the measured branching ratio, 
B(B° -> = (26.8 ± 1.3) x 10" 4 , this implies 

B(B° -> £> + Zr) T = -> D + ,D ) T = (3.64i°j 2 ) x 10~ 4 . (37) 

This result is in agreement, well within la, with the branching ratios measured by 
BaBar and Belle for B + — > D + D°. It deviates from the BaBar and Belle measurements 
of B(B° — ► D + D~) by 1.0a and 2.8a, respectively. (See Table HI) As we mentioned, 
deviations at this level are expected due to a term V* b V cd p ct and l/N c corrections. 

6.2.3 The amplitude p u and the smaller terms, pa u , e u , a u 

The CKM factor V* b V ud in AS = decay amplitudes multiplies four hadronic terms, 
p u ,pa u ,e u and a u . The QCD penguin amplitude p u = p ut , involving t and u quarks in 
the loop, is expected to have a magnitude comparable to that of p ct multiplying V* b V cd . 
Since we anticipate |V^,V^p rf | < 0.3|T| Eq. (1331) implies 

o.3 w: b v ud \ 



\v: b v udPu 




v; b v cd p ct \ 




\ V ub V ud\ 


\V*bVcdtc\ 




T + v; b v cdPct \ 


\v: b v cd \ 



< ■ ub "° = 0.2 . (38) 

" 1-0.3 \v: b v cd \ 1 ; 

This upper bound assumes a worst-case scenario of destructive interference between T 
and V* b V cd p c t in t c , as indicated by the factorization prediction (1341) which is larger than 
the corresponding experimental ratio (l3Tj) . 

The other three amplitudes, pa u , e u and a u , involving an interaction of a spectator 
quark, are expected to be suppressed by Aqcd/^b- For instance, the weak annihilation 
amplitude a u factorizes at leading order in 1/N C , 

a u = %Ci + ^)(0|Ml - 75Hfi + >( J D+D° M 0> + {l/N 2 c ) . (39) 

The |0) — > |F> + D°) matrix element of the 7 = 1 vector current is parameterized by one 
form factor, 

(D + (p D )D°(p D )\u 1 ,d\0) = fl I=1) (q 2 )(pn-PD)» ■ (40) 
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Combining the two factors in Eq. (1391) . one finds that the leading term in a u is propor- 
tional to the isospin breaking D meson mass difference which is negligibly small, 

au = ^(Ci + ^-)fBf { + =1 \m 2 B )(m 2 D+ - m 2 D0 ) + 0(1/N 2 C ) . (41) 

This implies that a u is dominated by nonfactorizable contributions including initial state 
gluon emission. 

Since pa u is also suppressed by a QCD loop factor, and e u and a u are suppressed by 
requiring a popping of a heavy cc pair, we will assume 

\pCL u | , \C U 

\,K\ « \Pu\ ■ (42) 

Thus, in the subsequent analysis we will neglect these very small amplitudes. [For the 
same argument, the interaction of a spectator quark and cc popping, we have neglected 
in Section [5] the two EWP contributions, PEew{c) and PAew(c).] 



7 Reiterating B —> DD decays and a bound on r 

Neglecting the very small amplitudes in (f42|) . Eqs. (|T6|) (|T8j) and the last of Eqs. (|T9l) 
become: 



a. 


A(B° - 


> D+D~) 


= V&Vcd {t c + e c ) + V: b V ud p, 


b. 


A(B° - 


^ D°D°) 




c. 


A(B + - 


-> D + D°) 




d. 


A(B° - 


- did;) 


= v; b v cdec , 


e. 


A(B S - 


* d + d;) 


= v* b v cd t c + v* b v udPu , 



(43) 

r = 5 = arg f ^) A 3 = . (44) 

\V: b V cd \ \t c + e c \ "I UeJ' V ; 

Using the two upper limits (|28l) and (|38il we find 

r < 0.3 . (45) 

We consider this a conservative upper bound, as it allows for a worst-case scenario of 
two destructive interference terms in the denominator of r and for a large enhancement 
of the penguin amplitude in its numerator. 

We will now discuss the rate and asymmetry measurements in B — > DD in light of 
the expressions ( 1431) . which are expected to hold to a very good approximation. First, we 
note that CP asymmetries in B° — > D°D° and B° — > DfD~ vanish in this approximation 
because the amplitudes for these processes involve a single CKM factor V* b V cd . [This 
justifies the discussion below Eq. (ITUI ).] The decay rates for these two processes, which 
are dominated by an amplitude e c , are equal in the SU(3) limit. A small decay rate 
difference between B° — ■> D°D° and B° — > D+D~ is expected due to two effects working 
in opposite directions, uu versus ss popping on the one hand and exclusive production 
of D°D° versus DfD~ on the other. 
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Second, as mentioned, neglecting e u and a u leads to AI — 1/2 in B — > DD implying 
the triangle amplitude relation (jS]). This relation was shown to be violated at 1.8a by 
branching ratios and asymmetries measured by Belle. The Standard Model AI = 3/2 
amplitude e u + a u is too small to account for an observable violation. 

The Belle asymmetries by themselves also show two unexpected features (see Table 
H]), which are not shared by the BaBar measurements: 

1. The above upper bound on r and the second Eq. (TS]) imply a theoretical upper 
limit on the direct asymmetry in B° — > D + D~ , 

\A CP (B° -> D + D-)\ < < 0.55 . (46) 

1 + r z 

The value measured by Belle is larger than this upper limit by 1.5cr. 

2. The asymmetries in B° — > D + D~ and B + — > D + D° are expected to be equal, up 
to second order corrections from an interference of the small amplitudes e c and p u . 
In contrast to this expectation, the Belle asymmetry in B° — > D + D~ is positive 
and large while the one in B + — > D + D° is very small. The difference between the 
two asymmetries involves a statistical significance of 3.6cr. The Standard Model 
interference of e c and p u and the amplitude pa u , which we neglected in P3h). are 
too small to account for this large difference between the two CP asymmetries. 



8 Determining 7 

The amplitude for B° — ■> D + D~ given in (fl6h ) or (HHa ) and a suitably chosen AS = 1 
amplitude in (1201) provide a sufficient number of observables for determining the weak 
phase 7 in the flavor SU(3) limit. This method has been proposed in Refs. [23] and [37]. 
Here we wish to recapitulate this method, showing that the determination of 7 in this 
way is very sensitive to uncertainties in SU(3) breaking. 

The two asymmetries S(D + D~) and Acp(D + D~) are given in Eqs. (T5]) in terms of 
the three parameters r, 5 and 7. We are assuming (3 = (21.5 ±1.0)°. A third equation for 
these parameters is provided by the ratio of CP-averaged decay rates for B° — > D + D~ 
and its U spin counterpart B s — > D+D~. Alternatively, one may use instead of the 
latter process the decay mode B° — > D^D~ or B + — > DfD°. In this case one would 
have to estimate the effect of the exchange amplitude e c contributing to B° — > D + D~ 
but not to the latter two AS = 1 decay processes. [See discussion above including the 
upper bound (|28|) .] 

Focusing our attention on the U spin pair B° — > D + D~ and B s — > D+D~ , we define 
an experimentally measured ratio of CP-averaged decay rates, 

p_ (Vcs\ 2 T(B ^D + D-) _ 1-\ 2 B(B°^D+D-) t s 

\Vj T(B S ^D+D-) A 2 B{B s ^D+D-)r ' 1 J 

where [1] t s /t = 0.939 ± 0.021 is the ratio of B s and B° lifetimes. Neglecting the 
second term in (T2"Uh) suppressed by the tiny CKM factor (|2"5"|) . and introducing an SU(3) 
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breaking parameter £ for the ratio of t c +e c amplitudes in B s — > DfD s and B° — ► D + D , 
one obtains: 

£ 2 ,R = 1 + 2rcos5cos7 + r 2 . (48) 

Thus, for a given value of £, the three observables S(D + D~), Acp{D + D~) and R in 
Eqs. (J2J) and (jlSl enable a determination of r, 5 and 7 up to discrete ambiguities. 

In order to obtain an analytic solution for 7, and to overcome two of its four discrete 
ambiguities, it is convenient to introduce another observable in B° — > D + D [38] (see 
also [37]), 

B ,|?£iW) (49) 

1 + \ ad+d-\ 

obeying with the two asymmetries, 

A 2 CP + S 2 + D 2 = 1 . (50) 
Expressing the new observable in terms of r, 5 and 7, 

D _ cos 2(3 + 2r cos S cos(2f3 + 7) + r 2 cos 2((3 + 7) 
1 + 2r cos 5 cos 7 + r 2 

we note that this quantity is positive for (3 = (21.5 ± 1.0)°, r < 0.3 [as required by fj45l) ]. 
and for arbitrary values of 5 and 7. This information on the sign of D, which remains 
undetermined by the two asymmetries using (1501 . avoids two discrete ambiguities in 7. 

Eqs. ([2]) and (j4"8l) and the positivity of -D imply the following equation for 7 in terms 
of A CP ,S and £ 2 i? [37]: 



+ v /l _ A 2 CP - £ 2 cos2(/? + 7) - Ssm2((3 + 7) - 1 _ _L_ 

cos2 7 - 1 £ 2 i? ' ^ ' 

The plus sign in front of the square root follows from the positivity of D. 

In order to demonstrate the high sensitivity of determining 7 to the value of the 
SU(3) breaking parameter £, we plot in Fig. 2(a) the dependence of 7 on £ for two sets 
of values for A C p, S and R: 

(1) BaBar central values, A C p = —0.11, S = —0.54, R = 0.52 (solid curve). 

(2) Central values of the averages of BaBar and Belle, A C p = 0.37, S = —0.75, R = 0.39 
(dashed curve). 

We do not use the Belle central values because they are non-physical, violating the 
inequality S* 2 + A 2 CP < 1. Also shown in the plot is a band describing the currently 
allowed la range for 7 [17], 7 = {67.6t%®)°. Fig. 2(b) shows values of r corresponding 
to solutions for 7 for case (1) (solid curve) and (2) (dashed curve). A solution obtained 
with too large values of r, r > 1.5 in case (1) and r > 0.6 in case (2), is not shown. 

The two curves of 7 for case (1) and (2) cross the allowed band for 7 at £ = 1.3 
and £ = 1.7, respectively. This difference in the values of the SU(3) breaking parameter 
follows from the large experimental errors in the B° — > D + D~ measurements. The 
slopes of the two curves at the above points are quite steep, implying in both cases a 
high theoretical sensitivity of the determined value of 7 to the assumed value of £. For 
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Figure 2: Behavior of solutions as a function of parameter £ describing SU(3) violation. 
Solid curves denote values obtained using BaBar central values as input; dashed curves 
based on BaBar-Belle average central values, (a) Weak phase 7. Horizontal dashed lines 
denote HFAG central and ±ler values of 7. (b) Parameter r = A u /A c describing ratio of 
amplitudes. Horizontal dot-dashed line denotes upper bound r < 0.3 described in text. 



instance, the dashed curve increases approximately linearly from 7 = 65° at £ = 1.7 to 
7 = 105° at £ = 1.6. Thus, assuming perfect measurements for the branching ratio and 
asymmetries in B° — > D + D~, an uncertainty of 10° in 7 would require knowing £ to 
better than 2%. A somewhat lower sensitivity has been noted in the second paper in 
Ref. [23]. 

The parameter £, representing SU(3) breaking in the ratio of t c + e c amplitudes in 
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B s — > DfD s and S — > D + D including a small calculable phase space effect, involves 
theoretical uncertainties from several sources: 

• As mentioned, t c include a dominant tree amplitude, for which the leading term in 
a 1/N C expansion factorizes. SU(3) breaking in the factorizable term is given by 
a product of a measured ratio of decay constants [36], /d s //d = 1-33 ± 0.07, and 
a ratio of form factors FB s ^£> s (m 2 Ds ) / ' Fb-*d{i , ^d) ■ The latter ratio is estimated at 
~ 1.05 in a chiral SU(3) perturbation expansion, at leading order in the heavy b 
and c quark masses m# [39]. However, a complete analysis of all 0(l/m#) terms 
shows that simultaneous violation of both chiral and heavy quark symmetries can 
be as large as 30% [40]. The form factor ratio can eventually be determined from 
B° — > D~£ + i> measured in e + e~ collisions on the T(4S) [35] and B s — > D~i + v 
accessible to the LHCb collaboration working at the LHC Thus, SU(3) breaking 
in the tree amplitude involves nonfactorizable 1/N C corrections and uncertainties 
in ratios of form factors and decay constants combining to a total of at least 30%. 

• A penguin contribution V* b V c dP c t in t c is not expected to factorize, which introduces 
an uncontrollable SU(3) breaking correction in this contribution at a level of 30%. 

• We have already discussed uncertainties in the magnitude of e c . 30% uncertainties 
due to SU(3) breaking corrections in this amplitude and in p ct mentioned above 
translate through (1251) and (I2"9"j) into two uncertainties in SU(3) breaking in t c + e c , 
each of which is at a level of 10%. 

In view of these combined uncertainties, a precision of 10% in £ is unachievable. This 
implies a very large intrinsic theoretical uncertainty of at least 50° in determining 7 as 
demonstrated in Fig. 2(a). 

The origin of the high sensitivity to SU(3) breaking can be traced back to the way 
flavor SU(3) symmetry is being applied here for a determination of 7. In this method 
an assumption of SU(3) symmetry is used to normalize the dominant amplitude in 
B — > D + D~ in terms of the measured amplitude for B s — > D+D~. In contrast, the 
error in 7 introduced by an uncertainty in SU(3) breaking is expected to be small in cases 
where a small penguin amplitude in AS = decays is normalized by a measurable SU(3) 
related AS = 1 decay amplitude. Two cases, where this has been demonstrated, are 
B° — > 7r + 7r~ [41], where the penguin amplitude is subdominant, and most prominently 
B° — > p + p~ [42] in which the penguin amplitude is very small. 

9 Conclusion 

We have studied the effect of a small b — > uud amplitude in B — > DD decays. While 
this amplitude violates AI = 1/2, we have argued that this violation is expected to be 
very small. Considering the Belle measurements, we pointed out a violation at a level of 
1.8cr of a AI = 1/2 amplitude relation, and an inconsistency at a level of 3.6a between 
CP asymmetries in B° — > D + D~ and B + — > D + D°. No such inconsistency has been 
observed by the BaBar collaboration. If these discrepancies persist they would have to 
be associated with a New Physics contribution to b — > uud. 
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Using conservative considerations within the CKM framework, we obtained a model- 
independent upper bound, r = A u /A c < 0.3, for the ratio of b — > uud and b — >■ ccd 
amplitudes in 5° — > D + D~. This implies an upper limit, \Acp\ < 0.55, for the direct 
CP asymmetry in this process. The Belle asymmetry is 1.5<r larger than this bound, 
while the BaBar measurement is well below this upper limit. 

U spin symmetry has been applied for obtaining relations between amplitudes and CP 
asymmetries for AS = decays of B and B s to pairs of charm-anticharm mesons and 
amplitudes and asymmetries in corresponding AS = 1 decays. We have shown that, 
while these relations are not useful for a precise determination of 7 in B° — > D + D~ 
in the presence of small uncertainties in SU(3) breaking, they may provide important 
information about small contributions to the latter process. 

For instance, the amplitude e c may account for a difference between the decay rates 
of B° — > D + D~ and B + — > D + D°, and for a small difference between CP asymmetries in 
these decays. This amplitude dominates the AS = decays, B° — ► D°D°, B° — ► D+D~ , 
and the AS = 1 decays B s — ► D + D~ , 5 S — > D°D°, and may be extracted from branching 
ratios of these processes. Using Eq. (127)) one would predict B(B° — ► D+Dj) = (4.0jr^;2) x 
10 -6 , whereas (125)) only implies £>(_B° — > D*D~) < 3 x 10~ 5 , close to the upper limit in 
Table I [15]. Similarly, using (J3D) B(B -> DfD) = (82.4 ± 6.5) x 10" 4 Eq. ([27]) would 
imply B(B S — > DD) = (0.9* x 10 -4 , while the conservative upper bound (125)) leads 
to a more modest prediction B(B S — > D/)) < 7 x 10 -4 . 

Thus, a sensitivity of 1 x 10~ 4 in B(B S — > D + D~) and B(S S — > D°D°), in comparison 
with — > D+D) = (82.4 ±6.5) x 10~ 4 , can be used for obtaining precise information 
on the ratio \e c /t c \ at a level of 0.1. This precision is considerably more powerful than the 
current upper bound ([281 obtained from the A 2 suppressed B(B° — > D°D°). The above 
sensitivity may be achieved at experiments carried by the LHCb collaboration working 
at the Large Hadron Collider. More precise measurement of the Cabibbo-favored decay 
branching ratio B(B S — ► D*D~) than currently available may soon be achieved at the 
Tevatron. This would lead to useful information on corrections to U spin symmetry 
relating this process and B° — > D + D~ . 
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A SU(3) operator analysis 

In this Appendix we consider relations between the graphical amplitudes defined in 
Sections 4 and 5 and a flavor SU(3) analysis for the operators appearing in the low 
energy effective Hamiltonian (jSJ). These operators have the following transformation 
properties under flavor SU(3). Current-current operators: 0% 2 ~ 3, C" 2 ~ 3,6,15, 
QCD penguin operators: O3-6 ~ 3, electroweak penguin operators: CV-10 ~ 3, 6, 15. 
An explicit SU(3) decomposition of the Hamiltonian can be found e.g. in [43]. 
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In B — ► decays, the initial and final states transform like 3 and 1,8, respec- 
tively, which implies that all these decay amplitudes can be expressed in terms of four 
SU(3) reduced matrix elements (1||3||3), (8||3||3), (8||6||3), (8|15|3). This agrees with 
the counting of independent amplitudes performed in Section 4 in terms of graphical 
amplitudes. 

The explicit expansion of all ten AS = and AS = 1 amplitudes in terms of reduced 
SU(3) matrix elements can be found, for example, in [44]. The AS = amplitudes are 
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(53) 



The corresponding AS = 1 amplitudes are given by the same transformation matrix, 
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(54) 



The subscript q = d, s on the reduced matrix elements is a reminder that they differ for 
AS = 0, 1 weak Hamiltonians through their dependence on CKM matrix elements. 

Comparing these expressions with the graphical expansions ( fT6l) and (!20il . we find 
the following relations between SU(3) reduced matrix elements and graphical amplitudes 
(for q = d,s): 



(11313), 

(81313), 

(81613), 
(8111513), 



1 1 

V*b v cq-^(tc + 3e c ) + V* b V uq -j=(p u + 3pa u + e u ) 



[8p u -e u + 3a u ) 



V* b V uq -(e u + a u ) 



(55) 



As mentioned in Section 5, electroweak penguin (EWP) contributions introduce four 
new graphical amplitudes. This agrees with the above counting of SU(3) reduced matrix 
elements. A complete expansion of EWP terms in the ten processes (1161) and ( 1201) in 
terms of graphical amplitudes defined in Section 5 is: 



Pew(B° - 
Pew(B° 



D + D~) 
* D°D°) 



v t tv td (^p£ w - 1 -pe 

2 2 
V t * b V td (--PA EW - -PA EW{c) ) 



2 1 N 

ew(c) + i^PAew — i:PAew(c)) 
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Pew(B + -> D + D° 

p EW (B° -> d+d; 



Pew{B s 
Pew{B s 
Pew(B + 
Pew{B s 
Pew(B° 



D + D~ 
DtD~ 



v: b v td {\p% w + ^pe 

V*V td (lpA RW - Ipa 



EW(c)) , 



\-EW(c)j 



VnYtdi^PEW 



-PE 



EW(c)j 



(56) 



2 G 1 2 1 

^t&^Mg-FW _ g-P-^£W(c) + -PA E w - -^PAew(c), 



v; b v ts {--PA l 

VttXtsi-PEW ~ 



3 

v^y to (|pA w - \pa 

v; b v ts {\p% w - \pe 



-PA EW ( C )) 
\PE, 



EW{c)) ■ 
EW{c)) 
EW{c)) ■ 



(57) 



Two of the EWP amplitudes can be related within a very good approximation to the 
amplitudes e u and a u appearing in Eqs. (IT5|) and (I2"0j) . Neglecting the EWP operators 
O7 8 which have negligibly small Wilson coefficients, SU(3) symmetry implies, 



PE 



EW(c) 



+ PA 



3 C 9 + C 



EW(c) 



10 



PE 



EW(c) 



-PA 



EW(c) 



2 d + C 2 

3 Cg — C\Q . 



(5? 



The proof of these relations is based on operator relations between 6 and 15 components 
of the EWP part of the effective Hamiltonian and corresponding components of the tree 
part (q = d, s) [43,45], 



n 



EWP 



(15) 



3 Cg + (7xo VttXtq 

"2 c!+c 2 v: b v uq 



74(15), «Wp(«) 



3 C 9 — C10 V^fcVfg 

2 d - C 2 KTfeKg 



?4(6) .(59) 



The first operator relation implies for q = s a, relation between AI = 1 EWP and tree 
amplitudes in B — > Kti decays [46]. 
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